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Abstract 

An  experimental  study  has  been  conducted  to  provide  insight 
into  heat  transfer  to  the  piston  of  a  light-duty  single-cylinder 
research  engine  under  Conventional  Diesel  (CDC), 
Homogeneous  Charge  Compression  Ignition  (HCCI),  and 
Reactivity  Controlled  Compression  Ignition  (RCCI)  combustion 
regimes.  Two  fast-response  surface  thermocouples  embedded 
in  the  piston  top  measured  transient  temperature.  A 
commercial  wireless  telemetry  system  was  used  to  transmit 
thermocouple  signals  from  the  moving  piston.  A  detailed 
comparison  was  made  between  the  different  combustion 
regimes  at  a  range  of  engine  speed  and  load  conditions.  The 
closed-cycle  integrated  and  peak  heat  transfer  rates  were 
found  to  be  lower  for  HCCI  and  RCCI  when  compared  to  CDC. 
Under  HCCI  operation,  the  peak  heat  transfer  rate  showed 
sensitivity  to  the  50%  burn  location. 

Introduction 

Increasingly  stringent  regulations  on  engine  emissions  and  fuel 
economy  have  motivated  research  to  improve  the  internal 
combustion  engine.  Currently  the  modern  direct-injected  diesel 
engine  is  the  most  efficient  engine  developed  for  automotive 
application.  The  direct  injection  of  diesel  fuel  relies  on  a 
mixing-controlled  combustion  event  with  a  diffusion  flame  that 
is  heterogeneous  in  nature  leading  to  fuel  rich  areas  that 
produce  particulate  matter  (PM)  and  stoichiometric  regions  that 
produce  nitric  oxides  (NOx)  [1 , 2].  NOx  and  PM  emissions 
have  been  shown  to  have  a  negative  impact  on  the 
environment  and  human  health.  Therefore,  government 
regulatory  agencies  limit  their  emission. 

Currently,  significant  research  has  been  focused  on  low 
temperature  combustion  (LTC)  regimes,  including  reactivity 
controlled  compression  ignition  (RCCI),  partially  premixed 
combustion  (PPC),  and  homogenous  charge  compression 
ignition  (HCCI).  Low  temperature  combustion  avoids  PM  and 
NOx  formation  by  keeping  the  local  equivalence  ratios  less 
than  unity  and  sufficiently  dilute  to  keep  the  temperature  below 
the  NOx  production  threshold  while  avoiding  PM  production.  In 
all  these  LTC  strategies  the  ignition  and  combustion  event  is 
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controlled  by  chemical  kinetics,  which  makes  combustion 
control  difficult  over  a  large  range  of  speeds  and  loads. 

HCCI  relies  on  a  completely  premixed  charge  of  air  and  fuel. 
The  ignition  process  is  kinetically  controlled  which  depends  on 
temperature,  pressure,  equivalence  ratio  and  fuel  properties  [3, 
4].  Researchers  have  reported  simultaneous  low  NOx  and  PM 
with  high  efficiencies.  However,  high  peak  pressure  rise  rates 
(PPRR)  limit  the  load  of  HCCI  due  to  practical  engine 
limitations.  Additionally,  HCCI  lacks  a  fast-response 
combustion  phasing  control,  such  as  spark  timing  or  injection 
timing,  when  a  single  fuel  is  used. 

RCCI  utilizes  in-cylinder  fuel  reactivity  stratification  by  utilizing 
a  second  fuel.  Typically,  a  low  reactivity  fuel  (such  as  gasoline) 
is  premixed  with  the  intake  charge  and  a  high  reactivity  fuel 
(such  as  diesel)  is  directly  injected.  The  spatial  variation  in  fuel 
reactivity  causes  a  sequential  auto-ignition,  reducing  the  PPRR 
relative  to  HCCI  combustion  [5-7].  RCCI  has  the  same  low 
NOx  and  PM  benefits  as  HCCI  but  has  the  added  benefit  of 
being  able  to  adjust  the  percent  of  directly  injected  fuel  on  the 
fly,  which  provides  control  of  the  combustion  phasing. 

The  in-cylinder  heat  transfer  has  been  studied  for  both 
convention  diesel  (CDC)  engines  and  engines  that  operate 
using  a  LTC  strategy  [8,  9].  The  heat  transfer  for  CDC  has 
been  studied  on  both  piston  surfaces  and  on  stationary 
locations  such  as  on  the  engine  head  and  cylinder  liners  [10, 
11].  More  recently,  measurements  have  been  made  on  the 
piston  surface  with  engines  operating  under  RCCI  and  CDC 
combustion  [12],  Measurements  made  on  the  piston  surface 
are  harder  to  make  than  those  on  fixed  locations.  Mechanical 
grasshopper  linkage  systems  have  been  used  to  connect 
piston-mounted  thermocouples  to  data  acquisition  systems  [8, 
9].  These  systems  have  a  limited  life  and  must  avoid  higher 
engine  speeds  due  to  vibration  concerns. 

Wireless  telemetry  systems  alleviate  many  of  the  problems 
encountered  with  a  mechanical  grasshopper  linkage  and  have 
been  used  to  acquire  piston  temperature  data  [8,  13,  14].  A 
wireless  telemetry  system  will  typically  outlast  the  useful  life  of 
most  fast-response  thermocouples.  The  bandwidth  of  many 
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wireless  systems  is  sufficient  to  make  crank-angle-resolved 
measurements. 

To  calculate  heat  flux  from  time-resolved  surface  temperature 
measurements,  a  one-dimensional  assumption  must  be  made. 
Periodic  boundary  conditions  are  applied  to  obtain  a  sustained 
solution.  Many  researchers  have  used  a  Fourier  series 
representation  to  recover  the  surface  heat  flux  [9,  1 5].  This 
method  relies  on  a  Fast  Fourier  Transform  (FFT)  and  is  many 
times  referred  to  as  such.  More  recently,  Hendricks  used  an 
inverse  solver  with  Tikhonov  regularization  to  calculate  surface 
heat  flux  [16].  The  inverse  solver  has  the  advantage  of 
additional  stability  when  analyzing  noisy  data  without  relying  on 
direct  filtering  of  the  data.  It  does,  however,  effectively  filter  the 
data,  see  ref.  [16]  for  details. 

This  experimental  study  seeks  to  understand  the  differences  in 
piston  heat  flux  for  HCCI,  RCCI,  and  CDC.  Operating  the 
engine  at  a  variety  of  conditions  allowed  a  detailed  analyses  to 
be  conducted.  The  results  presented  seek  to  add 
understanding  to  an  important  difference  between  LTC  and 
CDC  which  is  heat  transfer. 

Experimental  Setup 

Single-Cylinder  Engine  Laboratory 

Experiments  were  conducted  in  a  light-duty  single-cylinder 
engine  based  on  the  General  Motors/Fiat  JTD  1 ,9L  four- 
cylinder  diesel  engine.  A  stock  four  cylinder  head  was  installed 
atop  a  single  cylinder  block.  Valves  for  cylinder  #2  were 
activated  by  the  dual  overhead  cams.  The  valvetrain  for  the 
remaining  three  cylinders  was  deactivated  by  removing  the 
valve  rocker  arms.  Table  1  (adapted  from  Dempsey  et  al.  [5]) 
shows  the  specifications  for  the  engine. 


Table  1. Engine  specifications 


Enaine  Geometry 

Base  Engine 

GM  1 .9L  Diesel 

Compression  Ratio 

16.3 

Displacement  (Liters) 

0.477 

Stroke  (mm) 

90.4 

Bore  (mm) 

82 

Intake  Valve  Closing 

-132°aTDC 

Exhaust  Valve  Opening 

1 12°aTDC 

Swirl  Ratio 

1.5 -4.8 

Piston  Bowl  Type 

Stock  (Re-entrant) 

Port  Fuel  Injectors 

Included  Spray  Angle 

20° 

Injection  Pressure  (bar) 

2  to  10 

Rated  Flow  (cc/sec) 

<  10 

Bosch  Common  Rail  Injector 

Number  of  Holes 

7 

Hole  Diameter  (mm) 

0.14 

Included  Spray  Angle 

155° 

Injection  Pressure  (bar) 

250  to  1 000  bar 

The  research  engine  is  equipped  with  a  choked  orifice  air  flow 
metering  system  that  is  fed  with  building  air.  An  electrical 


heater  is  used  to  increase  the  temperature  of  the  air.  The  air 
system  includes  an  exhaust  gas  recirculation  loop  (EGR)  to 
provide  dilution  when  necessary.  The  EGR-air  mixture  enters  a 
heated  intake  surge  tank  to  dampen  wave  dynamics.  Two  port 
fuel  injectors  (PFI),  plumbed  with  separate  fuel  systems, 
deliver  up  to  two  different  fuels  to  create  the  premixed  charge. 
There  are  two  swirl-control  vanes  that  provide  an  adjustable 
swirl  ratio  between  1 .5  and  4.8.  The  engine  has  an  additional 
high  pressure  common  rail  fuel  injector  that  allows  for  direct 
injection  of  a  third  fuel.  After  leaving  the  engine  the  exhaust 
enters  a  surge  tank  to  dampen  wave  dynamics  and  enhance 
mixing.  Small  samples  of  the  exhaust  gas  are  acquired  by  an 
AVL  41 5S  smoke  meter  and  Horiba  gas  analyzers.  Figure  1 
shows  the  experimental  laboratory  schematically. 


4-cylinder  engine  head 
cylinder  #1,3,4  deactivated 


Figure  1.  Experimental  engine  layout 

The  engine  is  controlled  with  National  Instruments  LabView 
system.  The  fuel  injection  systems  are  controlled  by  Drivven 
PCI  controller  cards.  The  system  controls  and  records,  engine 
airflow,  emissions  data,  EGR  valve,  pressures,  low  speed 
temperatures,  and  other  critical  engine  parameters.  The  fuel 
mass  flow  rates  were  measured  separately  by  Coriolis  flow 
meters.  A  shaft  encoder  with  1440  pulses  per  revolution  was 
used  to  trigger  the  high  speed  data  acquisition  system.  In¬ 
cylinder  pressure  was  measured  with  a  Kistler  6125A  pressure 
transducer  in  series  with  a  Kistler  501 0  charge  amplifier. 

Cylinder  pressure  acquisition,  conditioning,  and  apparent  heat 
release  rate  (AHRR)  calculation  are  done  as  outlined  by  [5]. 
The  AHRR  uses  a  thermodynamic  first  law  analysis  with  an 
assumption  of  a  constant  ratio  of  specific  heats.  The  AHRR  is 
essentially  the  chemical  heat  release  minus  the  heat  transfer  to 
piston  and  combustion  chamber  walls. 

Wireless  Telemetry  System 

Temperature  data  were  transferred  from  the  piston  through  a 
IR  Telemetries  model  3051  microwave  transmitter  system.  The 
electromotive  force  (EMF)  generated  at  each  thermocouple 
was  converted  to  a  frequency  modulated  square  wave  over  the 
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range  of  10-50kHz.  A  voltage-controlled  oscillator  then 
converts  the  square  wave  into  a  microwave  signal.  The 
microwave  signal  is  broadcast  into  the  crankcase  where  four 
receiving  antennas  are  mounted.  The  signal  is  amplified, 
converted  back  to  a  modulated  square  wave,  downconverted, 
then  the  original  EMF  is  recorded.  The  signal  conversion 
process  induces  a  degree  of  filtering  but  the  advertised 
bandwidth  of  the  system  is  10kHz  [17].  Similar  systems  have 
been  used  by  other  researchers[12,  13]. 


a  plated  junction  formed  by  plating  the  surface  of  the  probe 
[18],  ora  sliver  junction,  formed  by  cold-weld  slivers  of  one 
material  onto  the  next.  Both  types  of  probes  were  utilized  in 
this  work.  The  thermocouples  were  J-type,  with  the  center  wire 
consisting  of  Constantan  and  the  outer  tube  of  Iron.  The 
thermocouples,  Medtherm  part  numbers  TCS-061-JU-0.25  and 
TCS-061  -JU-0.25-CR  were  claimed  to  have  a  time  constant  of 
one  microsecond  for  the  plated  junction  probes  and  a  ten 
microsecond  time  constant  for  the  sliver  junction  probes  [19]. 


The  microwave  transmitter  only  broadcasts  data  from  one 
channel  at  a  time;  it  cycles  though  the  seven  different  channels 
(+1  marker  channel).  Thermocouple  data  was  only  recorded  on 
two  channels  due  to  sensor  failure  (which  will  be  discussed 
later).  The  transmitter  was  configured  to  broadcast  each 
channel  for  0.5  seconds  independent  of  engine  speed.  At 
higher  engine  speed  more  engines  cycles  will  be  transmitted 
than  at  slower  engine  speeds.  The  data  streaming  off  the 
engine  is  collected  for  a  fixed  amount  of  time  and  the  data  are 
saved  to  a  local  computer. 

The  electronics  associated  with  the  wireless  telemetry  system 
are  powered  by  an  inductive  power  supply.  An  exciter  coil, 
mounted  on  the  underside  of  the  piston  liner,  is  powered  by  an 
inductive  power  supply.  When  the  piston  is  near  bottom  dead 
center  (BDC)  a  secondary  coil  that  is  mounted  on  the  piston 
skirt  interfaces  with  the  exciter  coil  allowing  power  to  transfer 
inductively.  Capacitors  and  other  electronics  store  enough 
energy  to  power  the  piston  though  the  entire  engine  cycle. 


Thermocouples 


Power 

Converter 

Primary  Coil 
(Engine 
Mounted) 


Inductive 

Power 

Supply 


Piston 


i-Transmitter 


Figure  2.  Wireless  Telemetry  System 

Thermocouple  Specifications 

Seven  fast-response  coaxial  thermocouples  were  installed  in 
the  piston  crown  to  measure  dynamic  surface  temperature. 
Thermocouple  data  was  only  recorded  on  two  channels  due  to 
sensor  failure.  Coaxial  thermocouples  have  been  used  in  a 
variety  of  applications.  The  sensor  can  survive  harsh 
environments  including  the  high  temperatures,  pressures  and 
accelerations  typically  encountered  in  an  engine. 
Thermocouples  with  both  plated  and  slivered  junctions  were 
installed  in  the  piston.  Surface  thermocouples  can  either  have 
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Table  2  and  Figure  3  list  and  depict  the  location  of  the 
thermocouples  mounted  in  the  piston.  All  the  thermocouples 
were  mounted  in  a  plane  perpendicular  to  the  wrist  pin.  Figure 
4  shows  a  top  view  of  the  piston. 

Upon  engine  startup  thermocouples  4,  5,  and  6  registered  no 
signal.  A  problem  with  the  thermocouples,  wiring  or  telemetry 
system  did  not  allow  any  data  to  be  collected  on  these 
channels.  Throughout  testing  thermocouples,  1  and  2  reported 
data  intermittingly.  It  is  believed  that  deposit  buildup  may  have 
caused  these  issues.  Thermocouples  3  and  7  provided 
repeatable  data  throughout  the  duration  of  testing. 

Table  2.  Thermocouple  Locations 


Number 

R  On) 

Z  (in) 

Angle 

(deg) 

Plating 

1 

0 

-0.217 

VERTICAL 

Chromium 

2 

0.3 

-0.425 

39.6 

Chromium 

3 

0.644 

-0.605 

VERTICAL 

Chromium 

4 

0.861 

-0.527 

40 

None 

5 

0.813 

-0.113 

60 

None 

6 

1.118 

-0.006 

VERTICAL 

None 

7 

1.31 

0 

VERTICAL 

Chromium 

Figure  3.  Piston  cross-section  showing  fast  response  surface 
thermocouples 
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Figure  4.  Piston  top  view  showing  fast  response  surface 
thermocouples 

Temperature  Data  Processing 
Data  Reduction 


The  temperature  data  recorded  by  the  wireless  telemetry 
system  included  spurious  spikes  that  needed  to  be  removed 
before  the  heat  flux  could  be  calculated.  The  spikes  in  the  data 
were  caused  by  the  wireless  telemetry  process.  When  the 
piston-mounted  transmitter  has  direct  line  of  sight  with  one  of 
the  receiving  antennas,  the  transfer  of  data  from  the  piston  to 
the  antennas  is  nearly  100%.  During  times  of  the  engine  cycle 
when  the  transmitter  does  not  have  good  visibility  with  the 
antennas  (i.e.  near  TDC),  the  transmission  efficiency  is 
degraded.  Figure  5  shows  the  raw  temperature  data  as 
recorded  by  the  vendor  software.  Useful  temperature  data  are 
buried  among  the  non-reasonable  points  that  include  negative 
values  and  values  beyond  the  melting  point  of  the  piston. 


Figure  5.  Temperature  data  including  spurious  spikes  before 
processing  on  a  sample  data  set 


Spurious  raw  data  were  eliminated  by  applying  the  generalized 
extreme  Studentized  deviate  (ESD)  test  at  each  crank  angle. 
The  ESD  test  is  used  to  detect  one  or  more  outliers  in  a  data 
set  that  approximately  follows  a  normal  distribution  [20].  This 
statistical  tool  is  preferred  over  other  tests  because  it  does  not 
require  that  the  number  of  outliers  be  predefined  before  the 
test  is  run  [21],  Typically  5%  of  all  data  was  removed  by  this 
procedure.  In  the  region  between  -20degATDC  and 
20degATDC,  approximately  25%  of  the  data  were  removed. 
Figure  6  shows  temperature  data  after  the  ESD  test  was 
applied. 


Figure  6.  Temperature  data  after  application  of  ESD  outlier  elimination 
on  a  sample  data  set 

After  the  ESD  test  was  applied  the  remaining  data  were 
averaged  at  each  crank  angle  degree.  Figure  7  shows  the 
resulting  temperature  plot 


Heat  Flux  and  T emperature 


X 

3 


03 

0 
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Figure  7.  Full  cycle  temperature  and  heat  flux  data 
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Instantaneous  Heat  Flux  Calculation 

The  surface  heat  flux  of  an  internal  combustion  engine  can  be 
represented  by  two  components:  the  steady-state  and  the 
transient  response.  The  steady-state  component  is  calculated 
from  a  time-averaged  surface  temperature  measurement  and  a 
back-side  thermocouple  at  a  known  depth  beneath  the  surface. 
The  transient  component  is  calculated  using  a  time-varying 
surface  temperature  measurement.  The  thermocouples  used  in 
this  study  only  included  surface  thermocouples  and  therefore 
only  the  transient  heat  flux  was  calculated. 

To  calculate  the  surface  heat  flux  a  one-dimensional 
assumption  must  be  made.  A  Fourier  analysis  is  utilized  to 
decompose  the  temperature  data  into  frequency  components. 

In  Equation  (1 )  The  Fourier  coefficients  A„  and  Bn  are 
calculated  by  applying  a  Fast  Fourier  Transform.  Tm  is  the 
mean  surface  temperature,  n  is  harmonic  number  and  w  is  the 
angular  velocity  of  the  engine.  The  FFT  method,  as  is  it 
commonly  known,  has  been  used  by  other  researchers  to 
calculate  surface  heat  flux  in  1C  engines  [9-11]. 

T(t)  =  Tm  +  ^J  An  cos(ncot)  +  Bn  sin(77frf)] 

Equation  (2)  is  the  periodic  transient  solution  for  the  surface 
heat  flux.  All  variables  are  known  from  the  Fourier  analysis  with 
the  exception  of  the  thermal  properties  of  the  material,  the 
conductivity,  k,  and  the  thermal  diffusivity,  a.  The  thermal 
properties  used  in  this  study  will  be  discussed  next. 

9,™„s  =  ^Za|^[(A,  +B„)c°s(ncot)-(All -Bn)sm(no)t)] 

To  calculate  heat  flux  thermal  properties  must  be  chosen.  The 
J-type  thermocouples  used  have  Constantan  as  the  inner  wire 
and  Iron  as  the  outer  tube.  The  piston  is  constructed  from  an 
aluminum  alloy.  But,  the  simplified  one-dimensional  analysis 
only  allows  for  a  single  set  of  properties  to  be  used.  Three- 
dimensional  finite  element  analysis  (FEA)  has  shown  that  the 
heat  flux  measurement  most  closely  matches  the  input  value 
when  the  material  properties  for  the  inner  wire  are  used  [22]. 

The  thermal  properties  for  Constantan  were  used  in  equation 
(2)  to  calculate  the  heat  flux.  For  each  case  the  time  average 
surface  temperature  was  used  to  calculate  the  conductivity  and 
thermal  diffusivity  based  on  measurements  made  by  Sasaki  et 
al.  [23].  The  thermal  properties  are  allowed  to  vary  slightly  to 
account  for  changes  in  operating  conditions  from  run  to  run, 
but  were  held  constant  as  temperature  varies  throughout  the 
engine  cycle.  Figure  7  shows  the  full-cycle  temperature  and 
heat  flux.  The  heat  flux  is  relatively  constant  at  a  negative 
value  through  the  intake  and  exhaust  part  of  the  cycle.  This 
negative  value  is  a  result  of  only  solving  for  the  dynamic  part  of 
the  heat  flux  solution.  During  the  compression  and  expansion 
part  of  the  cycle,  the  heat  flux  increases  due  to  the  increased 
gas  temperature  first  due  to  compression,  and  then  from  the 
high-temperature  combustion  products. 

To  have  an  understanding  of  the  total  energy  transfer,  the 
transient  heat  flux  was  integrated  over  the  closed  part  of  the 
cycle  from  intake  value  closing  (IVC)  to  exhaust  value  opening 


(EVO).  This  value  is  reported  hereafter  as  the  closed-cycle 
integrated  heat  flux. 

Results  and  Discussion 
Operating  Conditions 

The  engine  was  operated  at  four  different  load  and  speed 
points  (referred  to  hereafter  by  mode  number)  as  listed  in 
Table  3  using  RCCI,  HCCI,  and  CDC  combustion  strategies. 
The  combustion  phasing  (CA50),  gross  IMEP,  intake 
temperature,  and  intake  pressure  were  matched  for  the 
different  combustion  regimes.  Matching  load  and  CA50  was 
thought  to  give  the  best  comparison  between  the  different 
regimes  although  this  may  mean  a  particular  combustion 
regime  was  operating  at  a  non-optimal  point.  The  global  air- 
fuel  ratio  varied  between  the  different  modes  but  varied  less 
than  15%  across  the  different  combustion  regimes  for  a  given 
mode.  Mode  2  was  also  used  to  conduct  a  parametric  study  of 
CA50  and  swirl  effects.  The  CA50  and  swirl  conditions  run  for 
these  tests  are  listed  in  Table  4  and  Table  5,  respectively.  Swirl 
was  varied  by  adjusting  the  swirl-control  vanes  mounted  on  the 
intake  ports  of  the  engine. 

Table  6  lists  the  different  fuels  used  for  the  various  combustion 
regimes.  In  order  to  control  combustion  phasing  with  HCCI  two 
different  fuels  were  injected  using  two  port  fuel  injectors 
connected  to  different  fuel  systems.  Gasoline,  with  a  pump 
octane  number  (PON)  of  91 ,  was  used  in  one  PFI  and  a 
mixture  of  gasoline  and  n-heptane  was  used  in  the  other  PFI. 
N-heptane  is  an  extremely  reactive  fuel  that  functioned  to  lower 
the  octane  number  of  the  gasoline.  Individual  control  of  the  two 
PFIs  allowed  the  reactivity  (octane  number)  of  the  mixture  to 
be  changed  allowing  CA50  to  be  set.  CDC  was  operated  using 
F76,  which  is  a  military  specification  diesel  fuel  having  ignition 
quality  and  distillation  curves  similar  to  an  ultra-low  sulfur 
diesel  (ULSD).  During  RCCI  operation  F76  was  directly 
injected  and  91  PON  gasoline  was  port  injected. 


Table  3.  Engine  operating  conditions  at  different  loads  and  speeds 


Mode 

Mode 

Mode 

Mode 

1 

2 

3 

4 

Speed  (RPM) 

1490 

1900 

2300 

2300 

IMEPg  (bar) 

4.2 

5.7 

5.7 

8 

CA50  (degATDC) 

4 

5 

4.5 

8 

Swirl 

1.5 

1.5 

1.5 

1.5 

Intake  Temperature  (C) 

75 

50 

50 

35 

Intake  Pressure  (kPa) 

115 

130 

130 

188 

ERG  (%) 

0 

0 

0 

55 

During  CDC  operation  a  single  injection  was  used  for  all 
modes  with  the  exception  of  mode  1  which  included  a  small 
pilot  injection  to  decrease  the  ignition  delay  of  the  main 
injection.  The  length  and  timing  of  the  main  injection  was 
adjusted  in  order  match  the  desired  CA50  and  load.  The 
appendix  lists  the  injection  timing  used  for  each  case. 
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Table  4.  Engine  operating  conditions  for  CA50  sweep 


Mode  2 

Speed  (RPM) 

1900 

IMEPg  (bar) 

5.7 

CA50  (degATDC) 

3,  5,  7,  9 

Swirl 

1.5 

Intake  Temperature  (C) 

50 

Intake  Pressure  (kPa) 

130 

EGR  (%) 

0 

Table  5.  Engine  operating  conditions  for  swirl  sweep 


Mode  2 

Speed  (RPM) 

1900 

IMEPg  (bar) 

5.7 

CA50  (degATDC) 

5 

Swirl 

1.5,  2.5,  3.5,  4.8 

Intake  Temperature  (C) 

50 

Intake  Pressure  (kPa) 

130 

EGR  (%) 

0 

Table  6.  Fuels  used  for  different  combustion  regimes 


Regime 

Fuel 

HCCI 

91  PON  Gasoline  /  n-heptane 

RCCI 

F76/91PON  Gasoline 

CDC 

F76 

Combustion  Phasing  Effect 

A  CA50  sweep  was  conducted  on  mode  2  for  CDC,  HCCI  and 
RCCI.  Figure  8-10  show  the  heat  release  rate,  heat  flux  at 
location  3,  and  heat  flux  at  location  7,  for  the  three  combustion 
modes  respectively.  The  temperatures  listed  on  the  plots 
corresponds  to  the  mean  surface  temperature  for  that 
particular  case.  CA50  was  swept  from  3degATDC  to 
9degATDC  in  2  degree  increments. 

For  the  CDC  case,  as  CA50  was  advanced  the  peak  heat 
release  decreased  as  seen  in  Figure  8.  Injection  timings  closer 
to  TDC  have  shorter  ignition  delays  and  therefore  a  lower 
premixed  burn  fraction.  Both  location  3  and  7  showed 
advancement  in  peak  heat  flux  as  CA50  is  advanced.  Location 
3  peak  heat  flux  increases  with  advanced  timings  while 
location  7  does  not  show  a  trend.  Location  7  is  in  the  squish 
region  away  from  the  spray  plumes,  and  was  less  affected  by 
the  change  in  injection  timing. 
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Figure  8.  heat  flux  and  heat  release  for  CDC  CA50  sweep 

The  HCCI  data,  Figure  9,  are  seen  to  be  sensitive  to 
combustion  phasing.  As  CA50  is  advanced  the  peak  heat 
release  rate  increases  and  the  burn  duration  shortens.  The 
shorter  combustion  duration  corresponds  to  an  increase  in 
PPRR.  The  squish  land  data,  location  7,  becomes  similar  to 
the  CDC  data  at  the  advanced  combustion  phasing. 
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Figure  9.  heat  flux  and  heat  release  for  FHCCI  CA50  sweep 

The  RCCI  data,  however,  do  not  exhibit  the  same  trend  as  the 
HCCI  data  when  CA50  is  changed.  As  seen  in  Figure  10,  for 
RCCI  combustion  the  heat  release  does  not  vary  much  as 
CA50  is  advanced.  The  peak  heat  release  rate  increases 
slightly  but  the  general  shape  of  the  heat  release  is  unaffected 
by  combustion  phasing.  The  heat  flux  at  locations  3  and  7  also 
remains  relatively  constant  over  the  investigated  range  of 
CA50. 
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Figure  1 0.  heat  flux  and  heat  release  for  RCCI  CA50  sweep 

Figure  1 1  summarizes  the  effect  of  combustion  phasing  on 
CDC,  HCCI,  and  RCCI  by  showing  the  peak  and  closed-cycle 
integrated  heat  fluxes.  For  all  cases,  the  peak  and  closed-cycle 
integrated  heat  flux  for  CDC  is  larger,  by  approximately  a  factor 
of  two,  than  for  HCCI  and  RCCI  combustion.  The  HCCI  peak 
and  closed-cycle  integrated  heat  flux  values  consistently  lay 
between  the  RCCI  and  CDC  data.  Location  7  peak  heat  flux  for 
HCCI  increases  as  CA50  is  advanced.  For  the  most  advanced 
case  the  peak  heat  flux  of  HCCI  is  similar  to  CDC.  Location  3 
in  the  bowl  region  does  not  have  as  strong  of  a  trend  for  HCCI 
as  location  7. 


constant  at  5degATDC  for  each  case,  and  the  IMEPg  was 
maintained  at  5.7bar.  The  appendix  shows  the  injection  timings 
for  CDC  and  RCCI. 


Figure  12  shows  that  the  heat  release  rate  for  CDC  is  not 
strongly  affected  by  the  swirl  level  when  the  CA50  is  controlled. 
The  peak  heat  flux  at  locations  3  and  7  do  not  show  a  strong 
trend.  If  one  looks  at  the  motoring  part  of  cycle  (before  0  deg), 
however,  both  locations  show  an  increase  in  heat  flux  with 
increasing  swirl,  but  the  effect  is  much  more  pronounced  at 
location  7.  Because  the  injection  timing  was  changed  to 
maintain  combustion  phasing,  it  is  difficult  to  infer  how  the  in¬ 
cylinder  fuel  distribution  is  changing  with  the  swirl  level,  e.g., 
where  the  spray  plumes  impact  the  piston.  The  mean  surface 
temperature  for  both  locations  is  seen  to  monotonically 
increase  with  increasing  swirl. 
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Figure  1 2.  heat  flux  and  heat  release  for  CDC  swirl  sweep 

The  effect  of  swirl  on  HCCI  combustion  is  seen  to  be  similar  to 
CDC.  The  heat  flux  at  location  7  shows  a  similar  increase  in 
heat  flux  during  compression  prior  to  TDC  as  was  seen  in  the 
CDC  data.  However,  at  location  3  the  effect  prior  to  TDC  is  not 
significant.  There  is  no  discernible  trend  seen  in  the  peak  heat 
flux  with  respect  to  swirl  for  the  HCCI  data.  Again,  as  was  seen 
for  CDC,  as  swirl  was  increased  the  mean  surface  temperature 
also  increased  at  both  locations,  suggesting  an  increase  in  the 
steady  state  component  of  the  heat  transfer. 


Figure  1 1 .  Closed-cycle  integrated  and  peak  heat  flux  for  CA50  sweep 
for  location  3  and  7 


Swirl  Effect 

The  swirl  level  was  varied  from  1 .5  to  4.8  using  the  swirl 
control  vanes  mounted  on  the  intake  manifold  of  the  engine. 
Figure  12-14  show  the  effect  of  swirl  on  CDC,  HCCI,  and  RCCI 
combustion,  respectively.  The  combustion  phasing  was  held 
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Figure  13.  heat  flux  and  heat  release  for  HCCI  swirl  sweep 

As  swirl  increased,  the  peak  AHRR  was  found  to  increase 
slightly  for  the  RCCI  combustion  data,  which  can  be  attributed 
to  the  decreased  gradients  in  reactivity  because  of  the  higher 
mixing  rate.  At  location  3,  the  heat  flux  profile  is  largely 
unaffected  by  the  swirl  level.  But,  the  mean  surface 
temperature  is  seen  to  increase  with  the  swirl  level,  suggesting 
a  higher  net  heat  flux  to  the  piston.  This  is  believed  to  be  a 
result  of  the  fact  that  the  heat  flux  reported  here  is  the  dynamic 
heat  flux.  Larger  differences  may  be  expected  if  a  backside 
temperature  was  measured  and  the  steady  state  component 
added  to  the  solution.  At  location  7  there  is  an  increase  in  peak 
and  closed-cycle  integrated  heat  flux  as  swirl  is  increased, 
which  is  also  accompanied  by  an  increase  in  mean  surface 
temperature. 


matched  for  each  mode  in  order  to  make  a  comparison.  For  all 
cases  and  thermocouple  locations,  CDC  had  a  larger  peak  and 
closed-cycle  integrated  heat  flux  than  HCCI  or  RCCI 
combustion.  For  most  cases,  RCCI  had  the  lowest  peak  and 
closed-cycle  integrated  heat  flux. 

Mode  1  was  run  at  1490  RPM  and  4.2  bar  IMEPg  as  shown  in 
Figure  15.  The  peak  and  closed-cycle  integrated  heat  flux  for 
CDC  is  much  larger  than  for  HCCI  or  RCCI  combustion.  The 
HCCI  and  RCCI  heat  flux  profiles  are  very  similar.  The  lower 
engine  speed  allows  more  time  for  kinetically  controlled 
reactions  to  proceed.  The  AHRR  for  HCCI  and  RCCI  are 
similar  therefore  not  much  difference  is  expected  in  heat  flux  or 
mean  piston  surface  temperature. 
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Figure  15.  Mode  1  heat  flux  and  AFIRR 
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At  mode  2,  the  load  and  speed  were  increased  to  5.7  bar 
IMEPg  and  1900  RPM,  respectively.  In  Figure  16,  the  peak 
AHRR  for  mode  2  is  much  higher  for  HCCI  and  CDC  than  for 
RCCI.  Although  the  AHRR  is  similar  for  HCCI  and  CDC,  the 
heat  flux  is  very  different.  The  diffusion  flame  of  CDC,  with 
much  higher  local  temperatures,  causes  the  heat  flux  to  be 
much  greater  in  magnitude  and  duration. 


Figure  14.  heat  flux  and  heat  release  for  RCCI  swirl  sweep 


Combustion  Strategy  Effect 


To  better  understand  the  differences  in  piston  heat  transfer 
between  CDC,  HCCI,  and  RCCI  the  engine  was  run  over  a 
range  of  load  and  speed  points.  The  CA50  and  load  were 
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Figure  16.  Mode  2  heat  flux  and  AFIRFt 


Figure  18.  Mode  4  heat  flux  and  AFIRR 


Figure  17  shows  the  mode  3  (2300  RPM  and  5.7  bar  IMEPg) 
results;  the  HCCI  combustion  event  is  shorter  in  duration  than 
CDC  and  RCCI.  The  peak  heat  flux  at  location  3  for  HCCI  is 
greater  than  for  RCCI  but  less  than  CDC.  At  location  7  the  heat 
flux  profile  for  HCCI  begins  to  start  to  look  more  similar  to  CDC 
than  RCCI;  the  peak  heat  flux  is  nearly  the  same  in  magnitude. 


Mode  3 


- CDC 

—  HCCI 

—  RCCI 


-20  -10  0  10  20 
Crank  Angle  [deg] 


-20  -10  0  10  20  30  40 

Crank  Angle  [deg] 


-20  -10  0  10  20  30  40 

Crank  Angle  [deg] 


Figure  17.  Mode  3  heat  flux  and  AFIRR 


The  mode  4  (2300  RPM  and  8  bar  IMEPg)  results  are  seen  in 
Figure  1 8.  For  this  higher  load  point,  combustion  phasing  was 
retarded  to  a  CA50  of  8  degATDC  in  order  to  keep  PPRR 
under  control.  The  AHRR  for  HCCI  shows  a  shorter 
combustion  duration  and  increased  peak  AHRR  when 
compared  to  RCCI,  and  the  CDC  data  show  a  considerable 
diffusion-combustion  tail.  The  heat  flux  at  both  locations  is 
similar  for  HCCI  and  RCCI,  while  CDC  is  again  larger. 


Figure  19  and  Figure  20  summarize  the  closed-cycle 
integrated  heat  flux  for  the  four  different  modes  at  locations  3 
and  7,  respectively.  For  all  cases,  CDC  shows  a  significantly 
larger  heat  flux  than  HCCI  and  RCCI,  and  this  is  consistent 
with  the  mean  surface  temperatures  seen  in  Figure  15-18.  For 
all  cases  except  mode  1 ,  the  location  7  RCCI  case  had  a  lower 
closed-cycle  integrated  heat  flux  than  HCCI.  Mode  1  was  run 
at  a  lower  engine  speed  so  combustion  durations  for  HCCI  and 
RCCI  are  more  similar. 
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Figure  19.  Comparison  between  CDC,  FHCCI,  and  RCCI  closed-cycle 
integrated  heat  flux  at  location  3 
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Figure  20.  Comparison  between  CDC,  HCCI,  and  RCCI  closed-cycle 
integrated  heat  flux  at  location  7 

Summary/Conclusions 

1 .  RCCI  and  HCCI  had  a  lower  peak  and  closed-cycle 
integrated  heat  flux  than  CDC  for  all  cases.  The  LTC  strategies 
have  lower  in-cylinder  temperatures  than  CDC,  which  is 
considered  to  be  the  main  effect.  Locally  rich  areas  in  CDC 
also  lead  to  the  formation  of  soot,  which  may  also  impact  the 
heat  flux;  soot  radiation  is  not  expected  to  be  present  in  LTC. 
Modeling  by  Kokjohn  et  al.  discusses  the  origins  of  different 
heat  flux  for  CDC  and  LTC  and  also  shows  similar  differences 
[24]. 

2.  HCCI  peak  heat  flux  was  similar  to  RCCI  except  for  cases 
when  the  combustion  duration  for  HCCI  was  significantly 
shorter  than  for  RCCI,  causing  a  violent  combustion  event  with 
a  high  PPRR.  The  high  PPRR  events  may  have  an  effect  on 
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Definitions/Abbreviations 

AHRR 

Apparent  heat  release 
rate 

ATDC 

After  top  dead  center 

BDC 

Before  top  dead  center 

CA50 

Crank  angle  at  50%  to 
total  heat  release 

CDC 

Conventional  diesel 
combustion 

EGR 

Exhaust  gas  recirculation 

EMF 

Electromotive  force 

ESD 

Extreme  Studentized 
deviate 

FEA 

Finite  element  analysis 

FFT 

Fast  Fourier  Transform 

HCCI 

Flomogeneous  charge 
compression  ignition 

UNCLASSIFIED 


IMEPg 

Gross  indicated  mean 
effective  pressure 

LTC 

Low  temperature 
combustion 

NOx 

Mono-nitrogen  oxides 

PFI 

Port  fuel  injector 

PM 

Particulate  mater 

PON 

Pump  octane  number 

PPC 

Partially  premixed 
combustion 

PPRR 

Peak  pressure  rise  rate 

RCCI 

Reactivity  controlled 
compression  ignition 

RPM 

Revolutions  per  minute 

ULSD 

Ultra-low-sulfur  diesel 
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Appendix 

Injection  timings  for  CDC  and  RCCI 


CDC 

RCCI 

Injection 

Timing 

Injection 

Duration 

Injection 

Pressure 

Injection 

Timing 

Injection 

Duration 

Injection 

Pressure 

Combustion 
Phasing  Effect 

3degATDC 

-10 

0.531 

1000 

-56/ -36 

0.296/0.261 

500 

5degATDC 

-8 

0.53 

1000 

-59/ -39 

0.82/0.256 

500 

7degATDC 

-6 

0.532 

1000 

-63  /  -43 

0.565/0.252 

500 

9degATDC 

-4.5 

0.534 

1000 

-65  /  -45 

0.255/0.250 

500 

Swirl  Effect 

1.5  Swirl 

-8 

0.53 

1000 

-59/ -39 

0.281/0.256 

500 

2.6  Swirl 

-8 

0.531 

1000 

-58/ -38 

0.270/0.266 

500 

3.5  Swirl 

-8 

0.53 

1000 

-59/ -39 

0.265/0.266 

500 

4.8  Swirl 

-8 

0.53 

1000 

-59/ -39 

0.265/0.266 

500 

Combustion 

Strategy 

Effect 

Mode  1 

-26  /  -6 

.0.25/0.453 

1000 

-44 

0.286 

500 

Mode  2 

-8 

0.53 

1000 

-59/ -39 

0.810/0.256 

500 

Mode  3 

-11 

0.531 

1000 

-57/ -37 

0.286/0.261 

500 

Mode  4 

-10 

0.607 

1000 

-60/ -40 

0.278/0.250 

500 

*Note  if  two  values  are  listed  a  dual  injection  stagey  was  used 
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